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Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
units  as  follows: 


Multiply 

By 

To  Obtain 

acres 

4046.856 

square  meters 

feet 

0.3048 

meters 

knots  (international) 

0.514444 

meters  per  second 

miles  (U.S.  nautical) 

1.852 

kilometers 

1  Introduction 


Maalaea  Harbor  is  located  on  the  southwest  coast  of  the  Island  of  Maui, 
Hawaii,  approximately  8  miles*  south  of  Kahului  (Figure  1).  The  existing 
harbor  facility  consists  of  a  90-ft-wide,  12-ft-deep  channel  and  an  11.3-acre 
dredged  basin.  The  Harbor  is  surrounded  by  a  lOO-ft-long,  90-ft-wide  break¬ 
water  on  the  south  side  of  the  basin,  an  870-ft-iong  breakwater  on  the  east 
side,  and  a  3(X)-ft-long,  50-ft-wide  paved  wharf  on  the  north  side  of  the  basin 
(Figure  2). 

During  severe  wave  conditions,  Maalaea  Hartwr  experiences  harbor  surge 
and  navigation  difficulties  in  the  entrance  channel.  The  surge  results  from  the 
existing  configuration  and  alignment  of  the  harbor  entrance,  which  allow  direct 
wave  propagation  through  the  chaimel  opening.  Surge  problems  cause  naviga¬ 
tional  hazards  and  prevent  safe  berthing  in  some  portions  of  the  harbor 
(Lillycrop  et  al.  1993).  In  early  1980,  a  hydraulic  model  study  was  conducted 
at  the  Coastal  Engineering  Research  Center,  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  to  investigate  the  stability  of  the  various  break¬ 
water  cross  sections  considered  in  the  proposed  plans  of  improvement.  Details 
of  the  study  are  provided  in  Carver  and  Markle  (1981)  and  in  the  “General 
Design  Memorandum  for  Maalaea  Harbor  for  Light-Draft  Vessels”  (U.S.  Army 
Engineer  Division,  Pacific  Ocean  1980). 

While  the  improvements  are  aimed  at  the  navigation  and  economic  issues 
described  above,  it  is  important  that  the  proposed  designs  have  sound  technical 
bases  in  meeting  the  environmental  design  criteria  for  small  harbors 
(U.S.  Environmental  Protection  Agency  1985).  One  of  the  key  criteria  in 
measuring  the  physical  influence  of  water  on  the  aquatic  system  is  the  flushing 
time;  the  amount  of  time  that  it  takes  to  exchange  the  water  within  the  harbor 
with  the  receiving  water  body. 

In  estimating  flushing  time,  the  tidal  prism  technique  has  been  widely  used 
(Harleman  1966,  Dyer  1973,  Callaway  1981,  van  de  Kreeke  1983).  However, 
the  tidal  prism  method  has  two  severe  limitations  for  this  application:  (a)  it 
does  not  account  for  wind  effects,  which  are  important  factors  in  Maalaea 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  on 
page  viii. 
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Figure  2.  Existing  harbor  conditions 


Harbor,  and  (b)  it  assumes  that  the  receiving  water  body  is  well-flushed  within 
tidal  cycles.  As  pointed  out  by  Sanford,  Boicourt,  and  Rives  (1992),  the  tidal 
flushing  of  the  embayment  depends  not  only  on  the  embayment’s  channel  cur¬ 
rents,  but  also  on  the  coastal  currents  and  on  the  amount  of  mixing  that  occurs 
outside  the  embayment.  In  order  to  accurately  predict  the  flushing  time  in  a 
wind-dominant  water  body,  such  as  Maalaea  Harbor,  a  three-dimensional 
numerical  model  with  a  domain  covering  the  harbor-bay  system  was  invoked 
to  conduct  the  study. 

The  purposes  of  this  report  arc:  (a)  to  document  the  flushing  analysis 
methodology  using  numerical  model  techniques,  (b)  to  estimate  flushing  char¬ 
acteristics  in  Maalaea  Harbor  and  provide  necessary  information  for  assessing 
the  environmental  impact  resulting  from  proposed  improvement  plans.  The 
organization  of  this  report  is  as  follows:  Chapter  2  presents  the  numerical 
model  formulation,  including  the  governing  equations  and  the  solution  proce¬ 
dure.  Chapter  3  describes  model  calibration  and  validation  using  field  data 
collected  during  July  and  August,  1993.  Chapter  4  presents  the  scenario 
description,  the  hydrodynamic  results,  and  the  impact  of  the  proposed  plans  on 
flushing  characteristics.  Chapter  5  summarizes  conclusions  from  the  various 
parts  of  the  study. 
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2  Description  of  the 

Numericai  Hydrodynamic 
Modei 


The  numerical  hydrodynamic  model  used  in  this  study  is  CH3D-WES 
(Curvilinear  Hydrodynamics  in  T^ee  Dimensions- WES).  The  salient  feature 
of  CH3D-WES  is  its  capability  of  resolving  the  complex  geometries  in  the 
system  by  making  computations  on  the  boundary-fitted  curvilinear  planform 
grid.  The  model,  originally  developed  by  Sheng  (1986),  was  extensively  mod¬ 
ified  by  WES  to  enhance  the  model  capability  and  improve  the  performance  of 
the  model  code  (Johnson  et  al.  1991). 


Governing  Equations 


The  hydrodynamic  equations  used  in  CH3D  are  derived  from  the  classical 
Navier-Stokes  equations.  The  governing  partial  differential  equations  are 
based  on  the  following  assumptions: 


a.  The  hydrostatic  distribution  adequately  describes  the  vertical  distribu¬ 
tion  of  fluid  pressure. 

b.  The  Boussinesq  approximation  is  appropriate. 

c.  The  eddy  viscosity  approach  adequately  describes  turbulent  mixing  in 
the  flow. 

The  basic  equations  in  a  right-handed  Cartesian  coordinate  system  {x,y,z)  are: 

=  0  ( 

^  ^  ^ 
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(2) 


du  Buv  ^  duw  =  yv  _  J_ 

"5?  dx  dy  Bz~  Po 


a 

■5^ 


av  auv  av^  avw  _  ,  i  ap 

-5r-  p:^ 


a 


(3) 


(4) 


ar  ^  aur  avr  a^r 

"5?  ”5F  ^  az 


p  =  p  (r»s) 

where 

(«,v,w)  =  velocities  in  x,y,  and  z  directions,  respectively 
t  =  time 

/=  Coriolis  parameter  defined  as  2Q  sin  (j) 

Q  =  rotational  speed  of  the  earth 
0  =  latitude 


(5) 


(6) 
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p  =  density 
p  =  pressure 

AjjJCu  =  horizontal  turbulent  eddy  coefficients 
=  vertical  turbulent  eddy  coefficients 
g  =  gravitational  acceleration 
T  =  temperature 
S  =  salinity 

Equation  4  implies  that  vertical  accelerations  are  negligible  and  thus  the 
pressure  is  hydrostatic.  Various  forms  of  the  equation  of  state  can  be  specified 
for  Equation  7.  In  the  present  model,  the  foraiulation  given  below  is  used: 

p  =  P/(a  +  0.698P) 

where 

P  =  5890  +  387  -  0.3757^  +  3S 

a  =  1779.5  +  11.257  -  0.07457^  -  (3.8  +  0.017)5 


Nondimensionalization  of  Governing  Equations 

The  dimensionless  forms  of  the  governing  equations  are  used  to  facilitate 
relative  magnitude  comparisons  of  the  various  terms  in  the  equations.  The 
following  dimensionless  variables  are  used: 


(« *,  V  *,  w  *)  =  (u,  V,  wxjz^yu, 

(X*,  y*.  z*)  =(;c,  y, 

(t;,  Tp  =  (<, 
t*  =  tf 
C*  = 

P*  =  (P  -  Po)/(Pr  -  Po> 

T*  =  (J  -  TJKT^  -  7„) 

A*  =  AJA^ 

Kh  =  ^ 

<  =  KJK 
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where 


(x^,  Ty )  =  wind  stress  in  x  and  y  directions 
C  =  water  surface  elevation 

These  definitions  yield  the  following  dimensionless  parameters  in  the  govern¬ 
ing  equations: 


Vertical  Ekman  Number: 

Lateral  Ekman  Number: 

Vertical  Prandtl  (Schmidt)  Number: 
Lateral  Prandtl  (Schmidt)  Number: 
Froude  Number: 

Rossby  Number: 

Densimetric  Froude  Number, 
where 


E,  =  A,,  /fZ^ 

Eh  ~  ^Hr  ! 

Er,  = 

E^'h  ~  ^Hr  I  ^Hr 

Eo  =  IfX^ 

Fro  -FJyfl 
e  =  (p,  -  /  p^ 


(10) 


5, ,  ,  p, ,  ,  Afjr ,  -A,, ,  Kffr ,  and  are  arbitrary  reference  val¬ 

ues  of  the  salinity,  temjrerature,  velocity,  density,  etc. 


External-Internal  Modes 


The  basic  equations  (Equations  1-4)  can  be  integrated  over  the  depth  to 
yield  a  set  of  vertically  integrated  equations  for  the  water  surface  C  and  unit 
flow  rates  U  and  V  in  the  a:  and  y  directions.  Using  the  dimensionless  vari¬ 
ables  (asterisks  have  been  dropped)  and  the  parameters  previously  defined,  the 
vertically  integrated  equations  constituting  the  external  mode  are: 


IT 


(du 

■57 


+ 


dv] 


=  0 


(11) 
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As  will  be  discussed  later,  the  major  purpose  of  the  external  mode  is  to  pro¬ 
vide  the  updated  water-surface  field  for  input  to  the  internal  mode  equations. 

The  internal  mode  equations  from  which  the  3-D  velocity,  salinity,  and 
temperature  fields  are  computed  are; 
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dhS  ^  d  f  J.  95" 
dt  "Sz 


In  these  equations  h  is  the  thickness  of  an  internal  layer,  k+ll2  and  k-112  rep¬ 
resent  the  top  and  bottom,  respectively,  of  the  vertical  layer. 

Transformation  of  Governing  Equations 

The  CH3D  model  utilizes  a  boundary-fitted  or  generalized  curvilinear  plan- 
form  grid  which  can  be  made  to  confoim  to  flow  boundaries,  providing  a 
detailed  resolution  of  the  complex  horizontal  geometry  of  the  flow  system. 
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This  necessitates  the  transformation  of  the  governing  equations  into  boundary- 
fitted  coordinates  (^.Tj).  If  only  the  (x,y)  coordinates  are  transformed,  a  system 
of  equations  similar  to  those  solved  by  Johnson  (1980)  for  vertically  averaged 
flow  fields  is  obtained.  However,  in  the  CH3D  model  not  only  are  the  (j:,y) 
coordinates  transformed  into  the  curvilinear  system  but  the  velocity  is 
also  transformed  such  that  its  components  are  perpendicular  to  the  (^,ti)  coor¬ 
dinate  lines.  This  is  accomplished  by  employing  the  definitions  below  for  the 
components  of  the  Cartesian  velocity  (u,v)  in  terms  of  contravariant  compo¬ 
nents  u  and  v: 

U  =  XpU  +  x_v 

(19) 

V  =y^u  +  y^v 

along  with  the  following  expressions  for  replacing  Cartesian  derivatives 

f.  =  j 

(20) 
/,  -  i 

where /is  an  aibitrary  variable  and  J  is  the  Jacobian  of  the  coordinate  transfor¬ 
mation  defined  as 

J  = 

With  the  governing  equations  written  in  terms  of  the  contravariant  compo¬ 
nents  of  the  velocity,  boundary  conditions  can  be  prescribed  on  the  boundary- 
fitted  grid  in  the  same  maimer  as  on  a  Cartesian  grid  because  u  and  v  are 
perpendicular  to  the  curvilinear  ceU  faces  (e.g.,  at  a  land  boundary,  either  u  or 
v  is  set  to  zero). 

The  vertical  dimension  is  represented  through  the  use  of  what  is  commonly 
referred  to  as  a  sigma-stretched  grid.  The  vertical  depth  is  discretized  in  a 
fixed  number  of  layers,  each  layer  equal  in  thickness  to  a  fixed  percentage  of 
the  local  depth.  The  sigma-stretched  grid  is  then  transformed  to  a  fixed-space 
grid  where  the  computations  are  easily  performed. 

With  both  the  Cartesian  coordinates  and  Cartesian  velocity  transformed,  the 
following  boundary-fitted  equations  for  u,  v,  w,  S,  and  T,  to  be  solved  in  each 
vertical  layer,  are  obtained: 
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where 


<^11  =  4  *  yl 


Gi2  =  G21  = 

Similarly,  the  transformed  external  mode  equations  become: 


ac  o  aG  av  „ 
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+  Horizontal  Diffusion 

Equations  27-29  are  solved  first  to  yield  water-surface  elevations,  which  are 
then  used  to  evaluate  the  water-surface  slope  terms  in  the  internal  mode  equa¬ 
tions.  The  horizontal  diffusion  terms  are  quite  lengthy  and  thus  are  omitted  in 
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this  report.  Full  documentation  of  the  terms  is  presented  in  Johnson  et  al. 
(1991)  for  the  internal  mode  equations.  Similar  expressions  for  the  diffusion 
terms  in  the  vertically  averaged  equations  can  be  inferred  from  those  for  the 
internal  mode. 


Finite  Difference  Scheme 

Finite  differences  are  used  to  replace  derivatives  in  the  governing  equations, 
resulting  in  a  system  of  linear  algebraic  equations  to  be  solved  in  both  the 
external  and  internal  modes.  The  external  mode  solution  consists  ofjhe  sur^ 
face  displacement  and  vertically  integrated  contravariant  unit  flows  U  and  V. 
All  terms  in  the  transformed  vertically  averaged  continuity  equation  are  treated 
implicitly  whereas  only  the  water-surface  slope  terms  in  the  transformed  verti¬ 
cally  averaged  momentum  equations  are  treated  implicitly.  If  the  external 
mode  is  used  only  as  a  vertically  averaged  model,  the  bottom  friction  is  also 
treated  implicitly.  Those  terms  treated  implicitly  are  weighted  between  the 
new  and  old  time-steps.  The  resulting  finite  difference  equations  are  then 
factored  such  that  a  ^-sweep  followed  by  an  Tj-sweep  of  the  horizontal  grid 
yields  the  solution  at  the  new  time-step. 

Writing  Equations  11-13  as 


(30) 


du  ^  H  ^  ac 
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=  M 


(31) 


dV  ^  H 
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(32) 


(33) 
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A  typical  value  of  0  of  0.55  yields  stable  and  accurate  solutions. 


The  internal  mode  consists  of  computations  from  Equations  21-25  for  the 
three  velocity  components  (u,  v",  and  w),  salinity,  and  temperature.  The  only 
terms  treated  implicitly  are  the  vertical  diffusion  terms  in  all  equations  and  the 
bottom  friction  and  surface  slope  terms  in  the  momentum  equations.  Values  of 
the  water-surface  elevations  from  the  external  mode  are  used  to  evaluate  the 
surface  slope  terms  in  Equations  21  and  22.  As  a  result,  the  extremely  restric¬ 
tive  speed  of  a  free-surface  gravity  wave  is  removed  from  the  stability  criteria. 
Roache’s  second  upwind  differencing  scheme  is  used  to  represent  the  convec¬ 
tive  terms  in  the  momentum  equations,  whereas  a  spatially  third-order  scheme 
developed  by  Leonard  (1979)  (called  QUICKEST)  is  used  to  represent  the 
advective  terms  in  Equations  24  and  25  for  salinity  and  temperature,  respec¬ 
tively.  For  example,  if  the  velocity  on  the  right  face  of  a  computational  cell  is 
positive,  then  with  QUICKEST  the  value  of  the  salinity  used  to  compute  the 
flux  through  the  face  is 
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The  more  interested  reader  is  referred  to  the  paper  by  Leonard  (1979). 

It  should  be  noted  that  once  the  u  and  v  velocity  components  are  comput¬ 
ed,  they  are  slightly  adjusted  to  ensure  conservation  of  mass.  This  is  accomp¬ 
lished  byjorcing  the  sum  of  u  over  the  vertical  to  be  the  vertically  averaged 
velocity  UIH  and  the  sum  of  v  over  the  vertical  to  equal  V/H,  where  H  is  the 
total  water  depth. 
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r  - 
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Turbulence  Parameterization 

Vertical  turbulence  is  handled  by  using  the  concepts  of  eddy  viscosity  and 
diffusivity  to  represent  the  velocity  and  density  correlation  terms  that  arise 
from  time  averaging  of  the  governing  equations.  These  eddy  coefficients  are 
computed  from  mean  flow  characteristics  using  a  simplified  second-order  clo¬ 
sure  model  originally  developed  by  Donaldson  (1973).  The  closure  model  has 
been  further  developed  and  applied  to  various  types  of  flows  by  Lewellen 
(1977)  and  Sheng  (1982,  1986).  A  discussion  of  the  implementation  of  the 
turbulence  model  taken  from  Sheng  (1990)  follows.  For  more  details,  the 
interested  reader  should  refer  to  these  references  and  to  Johnson  et  al.  (1991). 

Assuming  local  equilibrium  of  turbulence,  i.e.,  no  time  evolution  or  spatial 
diffusion  of  the  second-order  correlations,  an  equation  relating  the  turbulent 
kinetic  energy  and  the  macro-scale  of  turbulence  to  the  mean  flow  shear  and 
stratification  (given  by  the  Richardson  number  Ri)  can  be  derived  as 

*  A[(.bs  *  3b  ^  lbh)Ri 
-  AbsO-  -  2b)]Q^+  6(5  +  3  +  4bs)Ri^ 

+  (65  -  A)(l  -  2b)Ri=  0 
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where 


Z?  =  0.125 


5=  1.8 


A  =  0.75 


A.\l{du/dzf  +  (3v73zf 


In  the  above  expression,  q  is  defined  as 


= 


v'v'  +  H-V/ 


and  A  is  the  macro-scale  of  turbulence.  The  quantities  u',  V,  and  w'  are  the 
turbulent  velocity  fluctuations  and  the  overbar  indicates  time  averaging. 

It  can  also  be  shown  that  the  following  relations  hold: 


A  1  ♦  “ 
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where 
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and 


(0 


OD 


bs 


(44) 


Thus,  after  the  velocity  shear  and  flow  stratification  are  determined,  q  can  be 
computed  from  Equations  38  and  39.  w'w'  is  then  determined  from 


1(1  -  2(0) 
2 


(45) 


Finally,  after  A  is  prescribed,  u'w'  and  can  be  computed  ftom  Equa¬ 
tions  40  and  41  and  the  vertical  eddy  coefficients  can  be  determined  from 
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q  A(1  -  (0) 
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In  addition  to  setting  A  =  0.65z  near  boundaries,  the  following  three  basic 
constraints  are  used  to  compute  A  at  a  vertical  position  r; 


£A 

dz 

A  < 

A  < 


<  0.65 


P 


Qcut 


(48) 

(49) 

(50) 


where  N  is  the  Bnmt-Vaisala  ftequency.  Equation  50  states  that  A  is  less  than 
a  fraction  of  the  spread  of  turbulence  as  measured  by  the  distance  between  the 
location  of  a  maximum  q^  to  where  ^  is  equal  to  25  percent  of  the  maximum. 
The  coefficient  is  on  the  order  of  0.15  to  0.25. 
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3  Implementation  of  the 
Hydrodynamic  Modei 


The  preceding  chapter  discussed  the  governing  equations  for  CH3D-WES 
and  the  finite  difference  scheme  used  to  solve  them.  In  this  chapter,  issues 
related  to  the  actual  implementation  of  the  model  will  be  addressed.  Three 
procedures  are  required  to  be  fulfilled  before  the  numerical  model  can  be 
enacted;  they  are;  (a)  the  generation  of  the  computational  grid  with  the  topo¬ 
graphic  information  assigned  to  it,  (b)  the  assigning  of  the  appropriate  initial 
and  boundary  conditions,  and  (c)  the  calibration  and  validation  of  the  model 
using  field-measured  data. 


Computational  Grid  and  Topography 

The  computational  grid  used  in  the  CH3D-WES  allows  boundary-fitted 
coordinates  to  fit  the  ship  channel  and  irregular  shoreline  with  accurate  and 
economical  grid  schematization.  With  the  aid  of  the  grid  generation  software 
(Jin  1993),  the  detailed  shoreline  configuration  including  the  breakwater  and 
the  coral  shoal  from  U.S.  Corps  of  Engineers  base  map  for  Maalaea  Harbor 
were  all  well  resolved,  as  shown  in  Figure  3.  Total  grid  cells  generated  in  the 
harbor-bay  system  were  9,720  with  81  in  the  east-west,  60  in  the  north-south 
and  2  layers  in  the  vertical  direction.  The  average  cell  size  is  50  ft,  with  the 
finest  resolution  measured  (25  ft)  being  mostly  concentrated  inside  the  harbor. 
The  computational  domain  was  bounded  by  the  shoreline  to  the  west  and  to 
the  north,  and  by  open  boundaries  in  the  east  and  the  south.  The  open  bound¬ 
ary  in  the  east  is  about  8(X)  ft  from  the  east  breakwater  and  the  open  boundary 
in  the  south  is  about  1,100  ft  from  the  harbor  entrance. 

The  topography  shown  in  Figure  4  was  read  from  the  1:800  scale  chart  and 
slightly  smoothed  through  the  discrete  digitization.  The  greatest  depth  is  25- 
30  ft  in  the  southeast  portion  of  the  computational  domain  and  is  shallowest  to 
the  east  of  the  east  breakwater  with  depths  of  2-4  ft.  The  depth  in  the  harbor 
basin  itself  is  generally  8-10  ft  except  for  a  ship  channel  located  from  the 
entrance  to  the  north  pier,  where  the  depth  is  14-16  ft. 
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Figure  3.  Curvilinear  grid  layout  for  Maalaea  Harbor  and  its  vicinity 


■ ;  4.42 
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Figure  4.  Maalaea  Harbor  and  Maalaea  Bay  topography 


A  vertically  stretched,  staggered  grid  scheme  is  used  internally  for  computa¬ 
tion.  Figure  5  illustrates  schematically  the  staggered  grid,  in  which  vertical 
velocity  is  shifted  one  half  an  interval  relative  to  the  density,  while  the  hori¬ 
zontal  velocity  grid  is  shifted  by  half  an  interval  horizontally  relative  to  den¬ 
sity.  The  vertically  stretched  grid  preserves  the  same  order  of  vertical 
resolution  for  the  shallow  and  deeper  part  of  the  water  body,  which  leads  to  a 
smooth  representation  of  the  topography. 
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Figure  5.  Vertically  stretched  and  staggered  numerical  grid 
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Initial  and  Boundary  Conditions 


When  initiating  a  run  of  CH3D-WES,  the  values  of  C,  u,  v,  w,  and  T 
are  set  to  zero  and,  values  of  salinity  and  temperature  are  read  from  input  files. 
Since  the  salinity  and  temperature  are  derived  from  prototype  measurements  at 
a  limited  number  of  locations,  inteipolation  and  smoothing  were  required  to 
produce  the  final  3-D  input  files.  Generally,  the  salinity  and  temperature  fields 
are  held  constant  for  the  first  few  days  of  a  simulatioa 

Boundary  conditions  at  the  fi'ee  surface  are 
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where 

C  =  surface  drag  coefficient 
W  =  wind  speed,  meters  per  second 
Boundary  conditions  at  the  bottom  are 
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where 


K  =  surface  heat  exchange  coefficient 
Tg  =  equilibrium  temperature 
Q  =  bottom  friction  coefficient 

Uj,Vj  =  values  of  the  horizontal  velocity  components  next  to  the  bottom 


Cj  is  given  by  Jin  (zj/z^)] 


(56) 


With  Zj  equal  to  one  half  the  bottom  layer  thickness, 

Zg  =  bottom  roughness  height 
k  =  von  Karman  constant 

The  surface  drag  coefficient  is  computed  according  to  Garratt  (1977)  as 
follows 

C  =  (0.75  +  0.067  WO  X  lO'^  (57) 


with  the  maximum  allowable  value  being  0.003.  The  surface  heat  exchange 
coefficient  K  and  the  equilibrium  temperature  are  computed  from  meteoro¬ 
logical  data  (wind  speed,  cloud  cover,  wet  and  dry  bulb  air  temperatures,  and 
relative  humidity)  as  discussed  by  Edinger,  Brady,  and  Geyer  (1974). 

Three  stations  were  used  to  collect  prototype  velocity  measurements  during 
the  July  27-August  4,  1993,  field  experiment  They  are:  stations  CMl  (inside 
the  harbor)  and  CM2  (at  the  entrance)  for  current  velocity,  and  gauge  TGI  (at 
the  entrance)  for  surface  height,  as  shown  in  Figure  6.  Figure  7  shows  tidal 
elevation  data  during  the  period  which  have  strong  semidiurnal  and  diurnal 
signals  ranging  from  2  to  2.5  ft.  Recorded  three  times  a  day,  the  wind  data 
obtained  at  0000,  0600,  and  1800  daily  are  the  actual  Maalaea  Harbor  wind 
data  telemetered  to  the  Kahului  Airport.  Figure  8  shows  wind  speed  and 
direction  data  during  the  experiment  period.  For  the  first  3  days,  the  wind  is 
from  the  northeast  with  speeds  ranging  between  5  and  20  knots.  For  the  rest 
of  the  experiment  period,  the  wind  is  predominantly  from  the  north  with 
speeds  ranging  from  2  and  25  knots.  During  the  period  July  28-30,  wind  data 
were  estimated  from  nearby  weather  stations  due  to  the  malfunction  of  the 
transmitter  at  the  Maalaea  weather  station. 
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Figure  6.  Current  meter,  tidal  gauge  locations  and  stations  for  model  output 
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Figure  7.  Tidal  gauge  data,  27  July  -  4  August  1993 


At  the  ocean  boundary,  the  water-surface  elevation  is  prescribed  along  with 
time-varying  vertical  distributions  of  salinity  and  temperature.  Specified  val¬ 
ues  of  salinity  and  temperature  are  employed  during  flood  flow,  whereas  dur¬ 
ing  the  ebb,  interior  values  are  advected  out  of  the  grid.  Normal  velocity, 
viscosity,  and  diffusivity  components  are  set  to  zero  along  sohd  boundaries.  In 
the  Maalaea  Harbor  application,  since  salinity  and  temperature  are  virtually 
unchanged  during  the  experiment  period,  a  constant  value  was  employed. 


Calibration  and  Validation  of  Hydrodynamic 
Model 


Numerical  modeling  of  hydrodynamics  and  transport  in  three  dimensions  is 
a  highly  complicated  task.  To  demonstrate  the  capability  of  the  model,  it  must 
undergo  calibration  and  validation.  Calibration  is  the  procedure  where  certain 
model  parameters  are  adjusted  to  maximize  the  agreement  between  model 
results  and  measured  field  data.  The  adjustable  parameters  in  this  model  are 
friction,  drag,  and  mixing  coefficient.  Once  the  calibration  procedure  is  com¬ 
pleted,  the  model  is  applied  without  further  adjustment  of  the  parameters  for 
validation.  Obtaining  good  comparisons  between  model  and  measured  data 

provides  confidence  that  the  model  can  simulate  the  hydrodynamic  condition 
in  the  study  area. 

The  Madaea  Harbor  numerical  model  calibration  was  conducted  over  a 
2.5-day  period  stai^g  at  0800  on  1  August  and  ending  at  1930  on  3  August, 
1993..  Measured  tidal  elevations  were  used  to  drive  the  open  boundary  condi¬ 
tion.  Wind  data  were  used  to  compute  the  surface  stress.  The  wind  direction 
showii  is  the  direction  from  which  the  wind  is  blowing.  The  angle  of  the 
direction  is  based  on  true  north  convention;  an  angle  equal  to  zero  defines  a 
wind  blowing  from  north  to  south,  and  angles  increase  in  a  clockwise  direc¬ 
tion.  Average  wind  speed  experienced  during  the  period  of  simulation  was 
4  m/sec  (8  knots),  with  winds  blowing  predominantly  from  the  northeast. 
Velocity  data  measured  at  the  entrance  and  inside  the  harbor  were  used  to 
calibrate  the  model. 

At  the  outset  of  the  calibration,  a  Manning’s  n,  the  friction  parameter,  equal 
to  0.035  was  used  for  the  entire  area.  However,  to  better  calibrate  against  the 
velocity  measurement  data,  a  Manrang’s  n  equal  to  0.040  was  assigned  to  the 
area  with  depth  less  than  .5  ft,  a  Manning’s  n  of  0.037  was  assigned  to  the  area 
having  depths  less  than  10  ft  and  greater  than  5  ft  while  n  equal  to  0.035  was 
used  for  the  rest  of  the  area.  Figures  9  and  10  show  the  final  comparisons  of 
channel  velocity  for  the  surface  and  bottom  layers  at  station  4.  Observation  1 
was  measured  every  2  min  by  single  depth  Endeco  current  meters  while  obser¬ 
vation  2  was  a  velocity  profiler  measurement  with  30  min  as  the  time  interval. 

A  good  agreement  between  modeled  and  measured  velocity  was  obtained  for 
both  data  sets.  Notable  exceptions  occurred  during  brief  periods,  such  as  at 
hours  15,  19,  and  22,  where  the  model  undeipredicts  the  surface  current 
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Figure  9.  Numerical  model  calibration;  model  and  data  comparison  for  surface  layer  at  station  4, 
1  August  -  3  August  1993 


STATION  4  (50,37)  BOTTOM 
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Figure  10.  Numerical  model  calibration:  model  and  data  comparison  for  bottom  layer  at  station  4, 
1  August  -  3  August  1993 


speed.  It  appears  that  during  periods  of  high  wind,  wind-generated  waves  can 
degrade  the  accuracy  of  the  Endeco  current  meter  readings. 

The  final  set  of  model  coefficients  thus  determined  is  as  follows: 


Manning’s  coefficient 
Horizontal  mixing  coefficient 
Base  vertical  mixing  coefficient 
Wind  drag  coefficient 


n  =  0.035-0.040 
Ayi=  1000  cm^/sec 
/ly/  =  10  cm^/sec 
C  =  (0.075-1-0.67  W)  X  10'^ 


Manning’s  coefficient  n  equal  to  0.035  was  assigned  for  the  entire  area,  except 
for  shallow  areas  having  depths  less  than  5  ft  but  greater  than  10  ft,  where  a 
Manning’s  n  equal  to  0.037  was  assigned,  and  for  areas  having  depths  less 
than  10  ft,  an  n  equal  to  0.040  was  assigned.  The  horizontal  mixing  coeffi¬ 
cient  of  1,(KX)  cmVsec  and  a  base  vertical  mixing  coefficient  of  10  cm^/sec 
were  used  throughout  the  simulation.  The  wind  drag  coefficient  is  calculated 
by  the  formula:  C  =  (0.75-1-0.67  W)  x  10'^,  where  W  is  measured  in  m/sec. 


Validation  of  the  Maalaea  Harbor  hydrodynamic  model  was  conducted  over 
a  5-day  period  starting  at  0730  on  27  July  and  ending  at  0730  on  1  August. 
Figures  7  and  8  show  the  overall  tide  and  wind  conditions,  which  are  similar 
to  the  calibration  period  except  that  tide  amplitude  is  somewhat  reduced  in 
magnitude.  The  data  used  for  validation  were  measured  velocity  data  at  sta¬ 
tion  4  and  station  2.  Current  measurements  at  station  4  were  taken  at  3  ft  and 
7  ft  above  the  bottom  over  a  water  depth  of  1 1  ft.  Model  and  data  comparison 
results  are  shown  in  Figures  1 1  and  12.  Predicted  velocity  components  match 
the  prototype  measurement  reasonably  well.  Comparing  surface  and  bottom 
velocities  reveals  that,  most  of  the  time,  the  flows  are  in  opposite  directions; 
i.e.,  the  surface  layer  flows  out  of  the  harbor  while  the  bottom  layer  flows  into 
the  harbor.  This  two-layer  flow  pattern  was  set  up  because  the  surface  layer 
moves  southward  under  influence  of  the  north  wind  while  the  bottom  layer 
forms  a  compensating  flow  in  the  opposite  direction.  At  station  2,  current 
measurements  were  taken  at  3  ft  above  the  bottom  over  a  water  depth  of  8  ft 
Figures  13  and  14  show  that  the  modeled  and  measured  north-south  (V)  com¬ 
ponent  velocities  agreed  reasonably  well.  For  the  east-west  (JJ)  component 
velocity,  there  was  a  lesser  degree  of  agreement  between  measured  and  pre¬ 
dicted  values.  For  instance,  between  day  1  and  day  1.5,  the  model  under¬ 
predicted  the  measured  data.  This  discrepancy  may  be  due  to  the  inaccuracy 
of  the  wind  data  and  the  surface  wave  effect  as  mentioned  earlier. 
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Figure  1 1 .  Numerical  model  validation:  model  and  data  comparison  for  surface  layer  at  station  4, 
27  July  -  1  August  1993 


STATION  4  (50,37)  BOTTOM 
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Figure  12.  Numerical  model  validation:  model  and  data  comparison  for  bottom  layer  at  station  4 
27  July  - 1  August  1993 


STATION  2  (32,56)  DEPTH  AVERAGE 
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Figure  13.  Numerical  model  validation:  model  and  data  comparison  for  depth-averaged  u  velocity  at 
station  2,  27  July  - 1  August  1993 


STATION  2  (32,56)  DEPTH  AVERAGE 
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Figure  14.  Numerical  model  validation:  model  and  data  comparison  for  depth-averaged  v  velocity  at 
station  2,  27  July  - 1  August  1993 


4  Scenario  Description  and 
Resuits 


In  the  preceding  chapter,  a  numerical  hydrodynamic  model  for  the  existing 
Maalaea  Harbor  was  developed,  calibrated,  and  validated.  In  this  chapter,  the 
model  will  be  used  further  as  a  project  design  tool  for  the  modified  plans 
proposed  in  the  Maalaea  Harbor.  The  model  output  thus  generated  can  be 
compared  to  that  from  the  existing  condition  to  assess  the  impact  resulting 
from  the  proposed  plans,  including  alternative  plan  2  (AP2),  alternative  plan  6 
(AP6),  and  modified  alternative  plan  2  (MAP2). 

Scenario  Description 

AP2  and  AP6,  shown  in  Figures  15  and  16,  are  the  focus  of  this  chapter  for 
hydrodynamic  and  flushing  impact  studies.  MAP2,  with  which  the  coral  shoal 
was  connected  by  a  concrete  causeway  to  the  mainland,  will  be  presented  sepa¬ 
rately  in  Appendix  A. 

Under  AP2,  changes  from  the  existing  harbor  are  as  follows: 

a.  A  620-ft-long  extension  to  the  existing  south  breakwater. 

b.  An  additional  400-ft-long  revetment  on  the  seaward  side  of  the  existing 
south  breakwater. 

c.  A  610-ft-long  entrance  channel  with  width  varying  from  150  to  180  ft 
and  depth  varying  from  12  to  15  ft. 

Under  AP6,  changes  are  as  follows: 

a.  A  600-ft-long  mole  structure  built  into  the  existing  harbor. 

b.  An  additional  berth  facility  built  on  the  inward  side  of  the  mole. 

c.  Removal  of  a  submerged  coral  reef  in  the  existing  haitor  to  provide 
more  berth  space  for  the  existing  south  breakwater. 
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Figure  15.  AP2 
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Figure  16.  AP6 


For  scenario  runs  with  the  proposed  plans,  most  of  the  input  files  stay  the 
same,  except  the  grid  and  topographic  files,  which  are  modified  to  reflect  the 
new  configurations.  The  forcing  functions  such  as  tide  and  wind,  and  calibra¬ 
tion  parameters  are  identical  to  those  used  in  the  existing  condition.  The 
period  for  the  simulation  begins  on  27  July  and  ends  on  1  August  1993.  For 
the  hydrodynamic  impact  study,  results  are  presented  in  terms  of  circulation 
pattern  (vector  plot)  and  current  velocity  (time  series  plot).  For  the  flushing 
impact  study,  results  are  presented  in  terms  of  contour  plot  and  time  series  for 
the  tracer  concentration.  The  four  stations  used  for  presenting  time  series  data 
are  shown  in  Figure  6. 

Analysis  of  Hydrodynamic  Impact  from  Proposed 
Plans 

Results  for  existing  condition 

In  general,  the  circulation  in  a  small  harbor  with  little  freshwater  input  is 
often  driven  by  the  tide.  However,  with  a  prolonged  steady  wind,  currents 
induced  by  wind  can  be  as  significant  as  those  from  astronomical  tides.  The 
circulation  in  Maalaea  Harbor  is  in  such  a  category.  While  its  tide  rises  up 
and  down  regularly  with  range  between  2  to  2.5  ft,  its  wind  blows  steadily 
from  the  north  and  northeast,  with  speed  varying  from  2  to  25  knots.  The 
prolonged  steady  winds  in  Maalaea  Harbor  are  the  result  of  trade  winds  tun¬ 
neled  by  the  two  mountain  ranges,  Puu  Lalua  (5,788  ft)  and  Puu  Kukui 
(10,023  ft),  just  north  of  the  harbor.^  The  funneling  effects  basically  rectify 
the  variable  trade  wind  into  a  steadier  north  and  northeast  wind  with  slightly 
magnified  speed,  as  evidenced  by  the  wind  record  in  this  region.  Figure  17 
shows  the  vector  plot  of  the  surface  layer  circulation  at  day  3  under  steady 
northeast  wind  at  flood  tide.  One  can  readily  observe  that  the  velocities  in  the 
bay  being  driven  by  the  northeast  winds  flow  in  the  southwestward  direction 
with  current  speed  varying  from  10-15  cm/sec.  A  small  portion  of  the  flow 
outside  the  east  breakwater  appears  to  have  a  reversal  due  to  the  interaction 
between  westward  flow  and  the  breakwater.  When  approaching  the  west 
coastline,  velocity  is  somewhat  diminished  and  the  direction  turned  toward  the 
south.  The  circulation  inside  the  harbor  is  such  that  a  clockwise  circulation 
around  the  coral  shoal  was  set  up  due  to  the  north  wind  pushing  the  water 
against  the  south  breakwater.  In  addition,  a  smaller,  counterclockwise  flow 
was  also  observed  in  the  west  end  of  the  harbor.  Overall,  the  velocity  inside 
the  harbor  is  2-5  cm/sec,  which  is  about  an  order  of  magnitude  smaller  than 
the  velocity  in  the  bay  (10-15  cm/sec). 

Figure  18  shows  time  series  plots  of  velocity  for  the  surface  and  bottom 
layers  at  the  entrance  of  the  harbor.  Of  note  is  the  two-layered  flow,  that  is 


*  Persoiud  Communication,  1994,  H.  Nakashima,  U.S.  Army  Engineer  Division,  Pacific  Ocean, 
Ft.  Shatter,  HI. 
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Figure  17.  Vector  plot  of  surface  layer  circulation  at  day  3  (0730,  30  July  1993)  for  the  existing 
condition 


V  VELOCITY  (cm/sec)  V  VELOCITY  (cm/sec) 


Figure  18.  Time  series  of  surface  and  bottom  channel  velocity  at  the 
entrance  of  the  harbor  for  the  existing  condition 
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the  opposite  velocity  direction  on  the  surface  and  at  the  bottom,  set  up  by  the 
near  steady  wind  from  the  north.  It  becomes  obvious  that  superimposed  on  the 
regular  astronomical  tide,  there  is  a  significant  vertical  shear,  set  up  by  a  pro¬ 
longed  north  wind.  This  steady,  wind-driven,  two-layer  circulation  is  what 
makes  the  dynamics  of  Maalaea  Harbor  circulation  unique  from  the  other  small 
harbors,  in  which  tidal  force  is  the  dominant  force. 


Results  for  AP2 

In  AP2,  a  breakwater  was  proposed  to  extend  from  the  existing  south 
breakwater  toward  the  southeast  direction  (Figure  15).  Figure  19  presents  the 
surface  layer  circulation  at  day  3  when  it  was  at  flood  tide  with  wind  from  the 
northeast  direction.  A  noticeable  effect  of  the  proposed  breakwater  on  the  bay 
circulation  was  seen  to  occur  at  the  south  side  of  the  new  harbor  entrance,  in 
which  the  proposed  breakwater  deflects  the  incoming  flow  from  its  original 
southwestward  direction  into  a  straight  southward  direction  with  accelerating 
flow  speed.  Associated  with  the  deflection,  the  original  southwestward  flow 
now  turns  northward  in  the  area  between  the  proposed  breakwater  and  the  west 
coastline.  Locally,  an  eddy  was  generated  around  the  tip  of  the  proposed 
breakwater,  exerting  a  circular  effect  on  the  flow  around  the  harbor  mouth. 

Inside  the  harbor,  the  circulation  pattern  is  similar  to  the  existing  condition; 
nevertheless,  the  velocity  magnitude  is  smaller  compared  to  the  existing  condi¬ 
tion  (by  10-15  percent).  The  exchange  rate  between  the  harbor  and  bay  water 
is  reduced  due  to  the  decrease  of  the  available  area  for  free  exchange.  Fig¬ 
ure  20  shows  time  series  of  the  surface  and  bottom  velocity  at  the  entrance  of 
the  harbor  for  AP2.  While  the  variation  of  velocity  pattern  is  similar  to  the 
existing  condition,  the  velocity  magnitude  is  smaller  by  10  percent.  Again,  the 
two-layered,  wind-driven  flow — ^that  is,  surface  layer  flows  outward  while 
bottom  layer  flows  inward — ^was  seen  to  persist  under  the  north  wind. 


Results  for  AP6 

In  AP6,  a  mole  was  proposed  to  be  installed  inside  the  existing  harbor  from 
the  south  breakwater  toward  the  nortteast  direction  (Figure  16).  Figure  21 
shows  surface  layer  circulation  at  day  3  when  it  was  at  flood  tide  with  wind 
from  the  northeast  direction.  Unlike  AP2,  the  proposed  mole  inside  the  harbor 
has  very  little  effect  on  the  outside  bay  water,  as  is  expected.  Inside  the  har¬ 
bor,  the  original  harbor  proper  was  divided  into  two  regions  by  the  mole. 
Inside  the  mole,  the  circulation  pattern  changes  only  slightly:  the  velocity 
magnitude  is  similar  to  the  existing  condition.  Outside  the  mole,  however,  the 
chaimel  velocity  at  the  harbor  entrance  was  reduced  dramatically  and  has  one 
of  the  smallest  values  among  all  conditions,  as  shown  in  Figure  22.  Compar¬ 
ing  time  series  for  all  three  conditions  (i.e.,  existing,  AP2,  and  AP6)  indicated 
that  velocity  for  AP6  was  oidy  30  to  40  percent  that  of  the  existing  condition 
(a  drop  of  60-70  percent),  as  indicated  in  Figures  23-24.  This  is  a  significant 
reduction  of  the  velocity  at  the  harbor  entrance  and  hence  can  have  profound 
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V  VELOCITY  (cm/sec)  V  VELOCITY  (cm/sec) 


Figure  20.  Time  series  of  surface  and  bottom  channel  velocity  at  the 
entrance  of  the  harbor  for  AP2 
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Figure  21 .  Vector  plot  of  surface  layer  circulation  at  day  3  (0730,  30  July  1993)  for  AP6 
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Figure  22.  Time  series  of  surface  and  bottom  channel  velocity  at  the 
entrance  of  the  harbor  for  AP6 
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Figure  23.  Time  series  of  channel  velocity  (surface  layer)  at  the  entrance  of  the  harbor  for  existing 
condition,  AP2  and  AP6 


STATION  4  (50,37)  BOTTOM 
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Figure  24.  Time  series  of  channel  velocity  (bottom  layer)  at  the  entrance  of  the  harbor  for  existing 
condition,  AP2  and  AP6 


influence  on  the  exchange  rate  between  the  harbor  and  the  bay.  The  reduction 
of  harbor  entrance  velocity  in  AP6  is  attributed  to  two  factors.  First,  the  cre¬ 
ation  of  a  constriction  between  the  mole  and  the  northern  berth,  which  acts  as 
a  control  section  limiting  the  velocity  through  it.  Second,  the  effect  of  the 
inverse  U-shaped  chaimel,  which  alters  the  pressure  distribution  around  the 
harbor,  and  hence  limits  the  exchange  of  harbor  water  and  bay  water. 


Analysis  of  Impact  on  Flushing  Characteristics 
from  Proposed  Plans 

Definition  of  flushing  time 

Flushing  of  a  semi-enclosed  water  body  can  be  defined  as  the  time  required 
for  a  conservative  tracer  C  to  decrease  to  36.8  percent  (1/e,  e  =  2.71828)  of 
its  initial  concentration  Cq.  This  is  the  widely  accepted  definition  for  e-folding 
time.  Given  a  fixed  amount  of  tracer  to  be  diluted,  the  longer  the  flushing 
time  required,  the  poorer  the  flushing  rate  is.  Conversely,  a  higher  flushing 
rate  will  result  in  a  shorter  time  to  exchange  the  prescribed  amount  of  water. 
For  adequate  flushing  of  a  basin,  Clark  (1983)  recommends  that  a  maximum 
time  of  2-4  days  should  be  safe  as  a  design  criterion  while  a  period  of  more 
than  10  days  should  be  considered  an  unacceptable  flushing  time. 


Method  of  calculating  flushing  time 

The  flushing  of  a  shallow,  semi-enclosed  basin  is  dynamically  related  to  the 
rate  of  exchange  between  the  basin  and  the  large  water  body  to  which  the 
basin  is  connected.  Primary  factors  that  influence  the  flushing  of  a  shallow 
basin  include:  astronomical  tide,  wind  force,  basin  geometry  and  topography, 
and  concentration  of  the  receiving  water  body. 

In  calculating  flushing  time,  the  tidal  prism  method  is  the  most  widely 
used.  However,  the  tidal  prism  method  is  limited  by  the  following  conditions: 
(a)  it  does  not  account  for  wind  effects;  (b)  it  assumes  that  the  receiving  water 
body  itself  is  well-flushed  within  tidal  cycles;  and  (c)  it  does  not  account  for 
complex  geometry  and  topography  in  most  cases.  In  the  present  approach,  all 
three  factors  are  considered.  Furthermore,  the  harbor  and  the  bay  are  con¬ 
structed  as  a  coupled  system  in  the  numerical  model  framework.  With  the 
three-dimensional  flow  field  available,  the  conservative  tracers  are  now  intro¬ 
duced  into  the  system  to  calculate  the  flushing  time  numerically.  The  equation 
used  to  calculate  the  concentration  is  an  advection-diffusion  equation,  similar 
to  the  one  used  for  calculating  salinity  (see  Chapter  2).  Both  are  without 
internal  source  and  sink  tenns. 
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Flushing  characteristics  In  Maalaea  Harbor 

Initially,  a  tracer  of  100  ppt  was  introduced  into  the  entire  harbor,  with  a 
0-ppt  tracer  concentration  in  the  bay.  The  boundary  condition  is  such  that  dur¬ 
ing  ebb,  interior  concentration  values  are  advected  out  of  the  grid,  while  during 
flood,  the  specified  return  value  will  be  applied.  As  time  passed,  the  tracer 
gradually  spread  over  a  large  area;  the  concentration  inside  the  harbor 
decreased  while  the  concentration  in  the  bay  increased  until  eventually  it 
reached  the  steady  state  when  the  two  are  equal. 

Figures  25-27  show  the  time  evolution  of  the  tracer  1  day  after  it  was  intro¬ 
duced  into  the  harbor.  Figure  25  shows  the  concentration  at  the  beginning 
time,  a  low  tide,  while  Figure  27  shows  the  concentration  at  the  end  time,  a 
high  tide.  Figure  26  shows  the  intemiediate  concentration.  While  the  contour 
usually  grows  in  an  irregular  shape  and  a  highly  complex  manner,  there  are 
consistent  patterns  recognizable  throughout  the  simulation  period.  First,  the 
highest  concentration  is  always  located  at  the  west  comer  of  the  harbor,  a 
region  at  the  far  side  of  the  harbor  entrance.  Second,  a  strong  concentration 
gradient  is  present  at  the  harbor  entrance,  separating  two  relatively  well-mixed 
regions  -  harbor  and  the  bay.  Third,  concentrations  at  the  northern  bank  of  the 
harbor  are  always  higher  than  those  at  the  southern  bank. 

Contour  plots  during  the  same  time  for  AP2  are  shown  in  Figures  28  and 
29  and  for  AP6  in  Figures  30  and  31.  Comparing  the  concentration  in  the  har¬ 
bor  for  aU  three  conditions  (existing,  AP2,  and  AP6)  reveals  that  the  highest 
concentration  occurred  under  AP6  (73.20  ppt),  followed  by  AP2  (47.70  ppt) 
and  the  existing  condition  (43.95  ppt).  The  fact  that  AP6  has  the  highest 
concentration  implies  that  it  has  the  poorest  flushing  of  all  three.  Further 
examination  of  the  contour  pattern  shows  that  AP6  also  has  the  highest  con¬ 
centration  gradient  near  the  constriction,  an  indication  that  the  narrow  section 
indeed  acts  as  a  control  section,  limiting  the  amount  of  water  allowed  to  be 
exchanged. 

Time  histories  of  the  tracer  were  recorded  at  stations  1,  2,  and  3  inside 
the  harbor.  Figure  32  shows  that  the  concentration  at  the  existing  condition 
takes  2.9  days  to  drop  from  its  initial  l(X)-ppt  to  37  ppt,  by  definition  the 
flushing  time.  Similar  calculations  show  that  flushing  time  is  3.3  days  for 
AP2  (a  14-percent  increase  compared  to  the  existing  condition)  and  6.3  days 
for  AP6  (a  117-percent  increase  compared  to  the  existing  condition).  Fig¬ 
ure  33  shows  that  at  station  2,  the  flushing  time  is  2.7  days  for  the  existing 
condition,  3.2  days  for  AP2  (an  18-percent  increase),  and  5.5  days  for  AP6  (a 
103-percent  increase).  Figure  34  shows  that  at  station  3,  the  flushing  time  is 
2.1  days  for  the  existing  condition,  2.3  days  for  AP2  (a  9.5-percent  increase), 
and  5.3  days  for  AP6  (a  152-percent  increase).  Table  1  summarizes  aU  the 
flushing  times  calculated  at  different  locations  under  various  conditions. 
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Figure  25.  Contour  plot  of  conservative  tracer  concentration  at  low  tide,  1  day  after  it  was  introduced  into  the 
existing  harbor  condition 


existing  harbor  condition 
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Figure  27.  Contour  plot  of  conservative  tracer  concentration  at  high  tide,  1  day  after  it  was  introduced  into  the 
existing  harbor  condition 
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Figure  28.  Contour  plot  of  conservative  tracer  concentration  at  low  tide,  1  day  after  it  was  introduced  into  the 
AP2  harbor  condition 
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Figure  29.  Contour  plot  of  conservative  tracer  concentration  at  high  tide,  1  day  after  it  was  introduced  into  the 
AP2  harbor  condition 
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Figure  30.  Contour  plot  of  conservative  tracer  concentration  at  low  tide,  1  day  after  it  was  introduced  into  the 
AP6  harbor  condition 
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Figure  31.  Contour  plot  of  conservative  tracer  concentration  at  high  tide,  1  day  after  it  was  introduced  into  the 
AP6  harbor  condition 


existing  stations 
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Figure  32.  Time  series  of  conservative  tracer  concentration  at  station  1 ,  covering  existing,  AP2,  and  AP6 
conditions 
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e  33.  Time  series  of  conservative  tracer  concentration  at  station  2,  covering  existing,  AP2,  and  AP6 
conditions 
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Figure  34.  Time  series  of  conservative  tracer  concentration  at  station  3,  covering  existing,  AP2,  and  AP6 
conditions 


Table  1 

Flushing  Time  Calculated  Under  Different  Conditions 

Station 

Existing  Condition 

AP2 

AP6 

Station  1 

2.9  days 

3.3  days 

6.3  days 

Station  2 

2.7  days 

3.2  days 

5.5  days 

Station  3 

2.1  days 

2.3  days 

5.3  days 

Average  flushing  time 
in  the  harbor 

2.6  days 

2.9  days 

5.7  days 

Longest  flushing  time 
in  the  harbor 

2.9  days 

3.3  days 

6.3  days 

From  the  above  analysis,  it  appears  that  AP2  imposes  a  moderate  impact  on 
the  flushing  time  characteristics  in  the  harbor,  an  overall  increase  of  13.7  per¬ 
cent  compared  to  the  existing  condition.  In  contrast,  AP6  imposes  a  much 
stronger  impact  on  the  flushing  time  characteristics,  an  overall  increase  up  to 
117  percent  when  compared  to  that  of  the  existing  condition. 
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5  Summary  and  Conclusion 


Based  on  the  results  presented  in  Chaptere  2  through  4,  conclusions  are 
drawn  with  regard  to  the  perfomance  of  CH3D-WES,  and  its  prediction  of 
hydrodynamic  and  flushing  properties  under  various  conditions. 

Computational  Grid 

A  curvilinear  boundary-fitted  coordinate  was  used  to  generate  the  computa¬ 
tional  grid,  which  provides  enhancement  for  fitting  the  irregular  shoreline  and 
discrete  digitizing  topography.  The  grid,  represented  by  81  x  60  cells  with  an 
average  grid  size  of  50  ft,  covers  Maalaea  Harbor  as  well  as  the  adjacent  bay. 
The  adoption  of  the  harbor-bay  system  eliminates  unjustified  assumptions  and 
allows  more  realistic  calculation  of  flushing  time  for  the  small  harbor. 


Performance  of  the  Numerical  Model 

Calibration  was  made  by  comparing  model  results  and  the  field-measured 
data  during  the  intensive  measurement  period.  The  best  fits  of  the  parameters 
to  the  observed  velocity  distribution  are  found  for: 

Manning’s  coefficient  n  =  0.035-0.040 

Horizontal  mixing  coefficient  =  1,000  cmVsec 

Base  vertical  mixing  coefficient  Ay  =  10  cm^/sec 

Wind  drag  coefficient  C=  (0.075+0.67  W)x  10"^ 

where  W  is  wind  speed  in  m/sec. 

During  the  validation,  the  model  calculations  are  conducted  without  further 
changes  and  adjustments  of  the  parameters.  The  model  results  compared 
favorably  with  the  field-measured  data,  reaffirming  the  reliability  of  CH3D- 
WES.  In  the  vertical  direction,  the  model  results  show  a  persistent  two-layered 
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circulation  driven  by  the  prevailing  north  wind,  a  unique  feature  observed  in 
the  Maalaea  Harbor. 


Hydrodynamic  Results 

The  existing  circulation  in  the  Maalaea  Harbor  and  its  vicinity  is  strongly 
influenced  by  the  prolonged  north  wind,  resulting  from  the  rectification  of 
trade  wind  by  two  mountain  ranges  north  of  the  harbor.  The  velocity  in  the 
bay  generally  flows  southwestward  with  a  speed  of  10-15  cm/sec.  TTie  flow 
inside  the  harbor  exhibits  a  large  clockwise  circulation  around  the  coral  shoal 
with  a  speed  of  2-5  cm/sec.  At  the  hart)or  mouth,  superposed  on  the  regular 
astronomical  tide,  is  a  two-layered  wind-driven  flow  set  up  by  the  prolonged 
north  wind. 

Under  AP2,  the  proposed  breakwater  has  an  impact  on  the  bay  circulation 
in  that  it  deflects  the  incoming  flow  fiom  southwestward  to  southward  direc¬ 
tion,  which  in  turn  drives  a  northward  compensation  flow  in  the  region  west  of 
the  proposed  breakwater.  Inside  the  harbor,  the  circulation  pattern  is  similar  to 
the  existing  condition;  however,  the  magnitude  of  velocity  is  10-15  percent 
smaller.  The  exchange  between  the  harbor  and  the  bay  is  somewhat  reduced 
because  the  total  available  area  for  exchange  at  the  mouth  is  decreased  with 
the  proposed  breakwater. 

Under  AP6,  the  proposed  mole  structure  has  very  little  effect  on  the  bay 
water  outside  the  harbor.  Inside  the  harbor,  the  original  harbor  proper  was 
divided  into  two  regions  by  the  mole.  Inside  the  mole,  the  circulation  pattern 
changes  only  slightly;  the  velocity  magnitude  is  similar  to  the  existing  condi¬ 
tion.  However,  the  velocity  at  the  harbor  mouth  (outside  the  mole)  was  only 
30-40  percent  that  of  the  existing  condition,  a  significant  reduction  in  terms  of 
velocity  magnitude.  The  reduction  of  the  velocity  is  attributed  to  two  factors. 
First,  the  creation  of  a  constriction  between  mole  and  the  northern  berth,  which 
acts  as  a  control  section  limiting  the  velocity  through  it.  Second,  the  effect  of 
the  inverse  U-shaped  channel,  which  alters  the  pressure  distribution  aroimd  the 
harbor,  and  hence  limits  the  exchange  of  harbor  water  with  bay  water. 


Analysis  of  Flushing  Time 

The  concentration  of  conservative  tracer  was  calculated  using  advection- 
diffusion  equations  in  the  numerical  framework.  Initially,  a  conservative  tracer 
of  1(X)  ppt  was  introduced  into  the  entire  harbor  and  a  0-ppt  concentration  was 
used  in  the  bay.  As  concentration  evolved  in  a  complex  manner  with  time,  the 
contour  was  found  to  follow  certain  patterns.  First,  the  highest  concentration 
is  always  located  at  the  west  comer  of  the  harbor.  Second,  the  concentration 
at  the  northern  bank  of  the  harbor  is  always  higher  than  at  the  southern  bank. 
Third,  a  strong  concentration  gradient  always  occurs  at  the  narrow  section  of 
the  harbor.  The  typical  calculation  after  1  day  of  simulation  shows  that  the 
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highest  concentration  inside  the  harbor  is  43.95  ppt  for  the  existing  condition, 
47.70  ppt  for  AP2,  and  73.30  ppt  for  AP6.  AP6  has  the  highest  concentration 
and  thus,  the  poorest  flushing  of  all.  Based  on  the  e-folding  fomiulation,  it 
was  shown  that  the  existing  harbor  has  a  flushing  time  of  2.9  days.  AP2  has  a 
3.3-day  flushing  time  (a  14-percent  increase)  and  AP6  has  a  6.3-day  flushing 
time  (a  117-percent  increase). 

Additional  scenario  runs  were  conducted  for  Modified  Alternative  Plan  2 
(see  Appendix  A  for  MAP2),  in  which  the  coral  shoal  was  connected  to  the 
northern  berth  with  a  concrete  causeway.  The  flushing  time  for  MAP2  is  4.4 
days,  slightly  higher  than  that  of  AP2  (3.3  days)  and  lower  than  AP6  (6.3 
days). 


Chapter  5  Summary  and  Conclusion 


65 


References 


Callaway,  R.  J.  (1981).  “Flushing  study  of  South  Beach  Marina,  Oregon.” 

J.  Waterways,  Port,  Coastal,  and  Ocean  Engin.  Div.,  American  Society  of 
Civil  Engineers,  107(2),  47-58. 

Carver,  R.  D.,  and  Maride,  D.  G.  (1981).  “Stability  of  mbble-mound  break¬ 
water  Maalaea  Harbor,  Maui,  Hawaii,  hydraulic  model  investigation,” 
Miscellaneous  Paper  HL-81-1,  U.S.  Army  Engineer  Waterways  Experiment 
Station,  Vicksburg,  MS. 

Qaik,  J.  R.  (1983).  Coastal  ecosystem  management:  A  technical  manual  for 
the  conservation  of  coastal  zone  resources.  John  Wiley  &  Sons,  New 
York. 

Donaldson,  C.  dup.  (1973).  “Atmospheric  turbulence  and  the  dispersal  of 
atmospheric  pollutants.”  Workshop  on  Micrometeorology.  D.A.  Haugen, 
ed.,  American  Meteorological  Society,  Science  Press,  Boston,  MA,  313-90. 

Dyer,  K.  R.  (1973).  Estuaries:  A  physical  introduction.  John  Wiley  and 
Sons,  New  York. 

Edinger,  J.  E.,  Brady,  D.  K.,  and  Geyer,  J.  C.  (1974).  “Heat  exchange  and 
transport  in  the  environment,”  Report  14,  EPRI  Publication 
No.  74-049-(X)-3,  prepared  for  Electric  Power  Research  Institute, 

Palo  Alto,  CA. 

Garratt,  J.  R.  (1977).  “Review  of  drag  coefficients  over  oceans  and  conti¬ 
nents,”  Monthly  Weather  Review  105,  915-29. 

Harleman,  D.  R.  F.  (1966).  “Pollution  in  estuaries.”  Estuary  and  coastline 
hydrodynamics.  A.  T.  Ippen,  ed.  McGraw-Hill,  Inc.,  New  York,  630-47. 

Jin,  K.  R.  (1993).  “Integration  of  radiation  boundary  condition  and  adaptive 
grid  technique  in  C!H3D  model  application,”  unpublished  paper  prepared  for 
U.S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 


66 


References 


Johnson,  B.  H.  (1980).  “VAHM — A  vertical  averaged  hydrodynamic  model 
using  boundary-fitted  coordinates,”  Miscellaneous  Paper  HL-80-3, 

U.S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 

Johnson,  B.  H.,  Kim,  K.  W.,  Heath,  R.  E.,  and  Butler,  H.  L.  (1991).  “User’s 
guide  for  the  Chesapeake  Bay  three-dimensional  numerical  hydrodynamic 
model,  ’  Techmcal  Report  HL-91-1,  U.S.  Army  Engineer  Waterways  Exper¬ 
iment  Station,  Vicksburg,  MS. 


Leonard,  B.  P.  (1979).  “A  stable  and  accurate  convective  modeling  procedure 
based  on  upstream  interpolation,”  Computer  Methods  in  Applied  Mechanics 
and  Engineering  19,  59-98. 

Lewellen,  W.  S.  (1977).  “Use  of  invariant  modeling.” //anrffeoojt  o/mr&M- 
lence.  W.  Frost,  ed..  Plenum  Press,  New  York,  Vol  1,  237-80. 

Lillycrop,  L.  S.,  Bratos,  S.  M.,  Thompson,  E.  F.,  and  Rivers,  P.  (1993). 

“Wave  response  of  proposed  improvements  to  the  small  boat  harbor  at 
Maalaea,  Maui,  Hawaii,”  Miscellaneous  Paper  CERC-93-4,  U.S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 

Sandford,  L.  P.,  Boicourt,  W.  C.,  and  Rives,  S.  R.  (1992).  “Model  for  estimat¬ 
ing  tidal  flushing  of  small  embayments,”  J.  Waterways,  Port,  Coastal,  and 
Ocean  Engin.  Div.,  American  Society  of  Ovil  Engineers,  118(6),  635-54. 

Sheng,  Y.  P.  (1982).  “Hydraulic  applications  of  a  second-order  closure  model 
of  turbulent  transport”  Applying  research  to  hydraulic  practice.  P.  Smith, 
ed.,  American  Society  of  Qvil  Engineers,  New  York,  106-19. 

- .  (1986).  “A  three-dimensional  mathematical  model  of  coastal, 

estuarine  and  lake  currents  using  boundary  fitted  grid”  Report  N.  585, 
A.R.A.P.  Group  of  Titan  Systems,  Princeton,  NJ. 

- .  (1990).  “A  simplified  second  order  closure  model  of  turbulent 

transport,”  unpublished  paper  prepared  for  U.S.  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  MS. 

U.S.  Army  Engineer  Division,  Pacific  Ocean.  (1980).  “General  Design  Memo¬ 
randum  and  Final  Environmental  Impact  Statement:  Maalaea  Harbor  for 
Light-Draft  Vessels,”  Maui,  Hawaii. 

U.S.  Environmental  Protect  Agency.  (1985).  Coastal  marinas  assessment 
handbook.  United  States  Environmental  Protection  Agency  (USEPA) 
Region  FV,  Atlanta,  Ga. 

van  de  Kreeke,  J.  (1983).  “Residence  time:  Application  to  small  boat  basins,” 
J.  Waterways,  Port,  Coastal  and  Ocean  Engin.  Div.,  American  Society  of 
Civil  Engineers,  109(4),  416-28. 


References 


Appendix  A 

An  Additional  Numerical  Model 
Run  for  Modified  Alternative 
Plan  2 


upon  completion  of  numerical  model  runs  for  alternative  plan  2  (AP2)  and 
alternative  plan  6  (AP6),  an  additional  scenario  run  for  a  modified  alternative 
plan  2  (hereafter  MAP2)  was  requested  and  was  conducted  in  the  Maalaea 
Harbor.  The  modification  was  focused  on  the  coral  shoal  region  located 
slightly  north  of  the  central  portion  of  the  harbor.  Previously,  AP2  treated  the 
coral  shoal  region  as  an  isolated  revetment,  completely  detached  from  the 
mainland.  In  the  MAP2,  however,  the  revetment  is  linked  to  the  northern 
berthing  area  through  a  concrete  causeway  (see  Figure  Al). 

Having  adopted  MAP2,  the  numerical  hydrodynamic  model  was  re-run 
using  the  initial  boundary  conditions  and  forcing  functions  identical  to  AP2  for 
a  period  of  5  days.  Model  results  are  analyzed  and  presented  in  terms  of  the 
velocity  field  and  spatial  and  temporal  variation  of  the  conservative  tracer. 
Figure  A2  shows  the  spatial  distribution  of  the  surface  velocity  field  at  day  3 
(0730,  30  July  1993)  when  the  wind  is  blowing  from  the  northeast  and  the  tide 
stage  is  at  high  tide.  It  is  obvious  that  the  new  concrete  causeway  in  MAP2 
has  a  blocking  effect  on  flows,  preventing  water  from  moving  through  the  old 
passway  between  the  revetment  and  the  northern  berthing.  Figures  A3  and  A4 
show  the  spatial  plot  of  conservative  tracers  1  day  after  they  were  introduced 
in  the  harbor  at  high  tide  and  low  tide,  respectively.  Figure  A5  shows  the  time 
variation  of  the  tracer  concentration  for  the  entire  5-day  simulation  period. 
Initially,  concentrations  for  all  three  stations  were  equally  high;  then  each 
station  decreased  gradually  according  to  the  exchange  rate  with  the  receiving 
water.  Based  on  the  e-folding  formula,  flushing  time  is  the  time  required  for 
concentration  to  reach  37  percent  of  its  original  value.  Table  Al  lists  flushing 
times  under  different  conditions  for  each  of  the  three  stations. 

Under  MAP2,  the  calculated  longest  flushing  time  for  Maalaea  Harbor  is 
4.4  days,  which  is  higher  than  that  of  AP2  (3.3  days)  and  lower  than  AP6 
(6.3  days);  however,  it  is  still  under  the  EPA’s  coastal  marina  criteria  (5  days). 
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Table  A1 

Flushing  Times  Under  Different  Conditions  Including  MAP2 

station 

Existing 

AP2 

AP6 

MAP2 

Station  1 

2.9 

3.3 

6.3 

4.4 

Station  2 

2.7 

3.2 

5.5 

4.1 

Station  3 

2.1 

2.3 

5.3 

3.8 

Averaged  flushing 
time 

2.6 

2.9 

5.7 

4.1 

Longest  flushing 
time 

2.9 

3.3 

6.3 

4.4 

As  for  individual  stations  in  the  harbor,  station  1  in  the  far  west  has  the 
longest  flushing  time  (4.4  days),  followed  by  station  2  in  the  north  (4.1  days), 
and  station  3  in  the  south  (3.8  days).  The  average  flushing  time  is  4.1  days 
forMAP2. 
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Figure  A1 .  Geometry  for  MAP2 
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Figure  A2.  Vector  plot  of  surface  layer  circulation  at  day  3  (0730,  30  July  1993)  for  MAP2 


Figure  A3.  Contour  plot  of  conservative  tracer  concentration  at  low  tide,  1  day  after  tracer  was  introduced 
into  the  harbor  under  MAP2 
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MODIflED  ALTERNATIVE  PLAN  2 _ STATION  1  MODIFIED  ALTERNATIVE  PLAN  2  STATION  2 


Appendix  A  An  Additional  Numerical  Model  Run  for  MAP2 


A7 


Figure  A5.  Time  series  of  conservative  tracer  concentration  under  MAP2  for  stations  1 ,  2,  and  3 
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